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The dark matter constraint imposed by the recent WMAP experiment on gravitational theories is 
analyzed. Using the non-linearity of the vacuum Einstein's equations, it is shown that when the slow 
motion condition is applied to the geodesic equations, the resulting nearly Newtonian gravitational 
field describes nearly flat velocity curves for rotating stars in the vicinity of thin disk galaxies. 

I. INTRODUCTION 

The discrepancy between the observed rotation velocity curves for stars in a spiral galaxy and the theoretical 
prediction from Newtonian mechanics is a long standing problem in modern astrophysics. The velocity measurements 
are based on the Tully-Fischer relation between the mass of the galaxy and the width of the 21-cm line of hydrogen 
emissions, suggesting that a larger galaxy mass would increase the rotation rate. Since these velocities are much 
smaller than the speed of light c, in principle they should be described by Newton's gravitational theory. However, 
as shown in the example of Figure 1 for the galaxy NGC3198, the theoretical Newtonian curve agrees with the 
experimental one only at the galaxy's nucleus For larger distances, the observed curve becomes almost horizontal, 
separating from the theoretical Newtonian curve which drops rapidly with r. Such pattern is observed in most spiral 
galaxies and galaxy clusters 0, • 



100 



< 1 1 1 1 1 I 

10 20 30 

radius/kpc 

FIG. 1: The Observed Rotation velocity curve (error bars) compared with the predicted Newtonian curve of the NGC3198. 

The most common explanation for this problem was originally proposed by Zwicky in 1933 Q. Accordingly, a 
certain quantity of dark matter, invisible with respect to the electromagnetic radiation spectrum, should be added 
to each galaxy. Such matter can in principle be composed of ordinary baryonic matter, like planets distributed in a 
spherical halo orbiting the galaxy itself, far away from the stars [E| . These have been observed with the help of the 
gravitational microlensing effects, but only in very small amounts, far beyond the required quantity to correct the 
velocity curves. In the cosmological scale, dark matter seems to be consistent with the standard FRW model but 
only recently the cosmic microwave radiation data analysis from the WMAP experiment indicated that most of the 
estimated 22% dark matter content of the universe must be of non-baryonic nature. More specifically, the analysis of 
the power spectrum indicates that a theory of gravity based essentially on the properties of baryonic matter would 
produce a lower third peak Therefore, either some exotic particles must be considered Q, or else an adequate 
gravitational theory should be devised. 

In principle that constraint does not exclude non-linear theories like general relativity. However, general relativity 
has such strong commitments with its Newtonian limit, that it makes it difficult to explain the rotation velocity curves. 
The usual argument goes as follows: The velocity curves for stars in a spiral galaxy derived from Newton's theory 
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agree with the observed curves only at the galaxy's core (FIG. 1.), precisely where the space-time curvature produced 
by Einstein's gravity would be more pronounced. Beyond that region, the gravitational field becomes sufficiently weak 
to be taken over by its Newtonian limit. Over the time, this has motivated research on many alternative gravitational 
theories, which we separate into two main categories: 

(i) Modifications of Newtonian Gravity 

According to this proposal, Newton's gravitational theory should be modified so as to correctly describe the velocity 
curves. The first thought is of course the post-Newtonian approximations of general relativity, regarded as corrections 
to Newtonian theory. However, a simple exercise shows that in a second order parametric post-Newtonian approxima- 
tion, the correction term in the velocity curves decay with 1 /r 3 , not improving substantially the velocity curves |10| . 
On the other hand, post-Newtonian cosmology (see e.g. ^lj) does not meet the WMAP power spectrum constraint. 

Other modifications of the Newtonian theory have been considered pi IT2IIT3I ]. Among these, MOND has received a 
substantial attention and it has been backed by a theory in which Poisson's equation for the Newtonian gravitational 
field is replaced by an equation like Q 

< V,nC^-)V(p >=4ttG P 
a 

where fJb[x) is a function to be adjusted to the specific type of galaxy and ao is a constant acceleration. For example, 
in a spherically symmetric distribution of matter, it is suggested that n(x) — xj\Jx 2 — 1, producing the following 
correction for Newtonian potential 

ip = \J aoGM In r 

This theory has shown good agreement with most known spiral galaxies, but there are reports suggesting that it may 
be constrained by galaxy clusters [TH . Finally, in spite of being essentially a local gravitational theory, its global 
effect on the composition of the total energy of the universe does not meet the power spectrum constraint p| • 

(ii) Modifications of General Relativity 

There are just too many ideas on how to modify general relativity to correct the velocity curves based on a variety 
of suppositions. Here we just list some of these: (1) Add a scalar field to Einstein's equation, in such a way that the 
scalar-tensor theory corrects the Newtonian limit ; (2) Modify the concept of time in general relativity, so that the 
Newtonian limit of the theory differs from the original Newton's' theory [I&ll9| : (3) Add a cosmological constant with 
the appropriate sign (depending on which side of the equation it is placed) |2Qj ; (4) Include higher order curvature 
terms in the gravitational variational principle as a means to increase the local gravitational pull on galaxies |2lj |: 
(5) Several brane-world models and variants have been considered, in the hope that the additional degree of freedom 
would explain the rotation curves, pi I2I I2I I2I IM I2H |2S | . 

The purpose of this paper is to show that when the slow motion condition v << c is applied to the geodesic 
equations only, then the self interacting vacuum gravitational field produced by a disk galaxy, contributes to a nearly 
Newtonian motion of a star in the galactic plane, with nearly flat velocity curves. 

This is justified first by the fact that the geodesic equations are derived from Einstein's equations, but in the 
Newtonian limit the equations of motion corresponds to a separate postulate 0. Therefore, when the gravitational 
field of a galaxy acting upon a free falling star is sufficiently weak, then the slow motion condition v << c applies 
to Einstein's equations and the only remaining option is the Newton's law of motion. On the other hand, admitting 
that the free falling star gets in a region where the gravitational field is beyond the Newtonian limit, then the 
condition v << c still applies to the geodesic equations, but not necessarily on Einstein's equations. In fact, the 
geodesic equations depend only linearly in the connection, while Einstein's equations depend quadratically in the 
same connection. Therefore, the effect of the condition v << c in the connection, becomes less restrictive on the 
geodesic equation than in Einstein's equations. The result is that in that region a nearly Newtonian gravity prevails. 

II. NEARLY NEWTONIAN GRAVITY 

Consider a slowly free falling star, u « c, in a region of the space-time, where the pull of the gravitational field 
on the particle is initially weak: 

9nv='rinu+5h liv , dh 2 ^ « Sh^, (1) 
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where 5h^ u is not parameterized by v/c. Under these conditions Newtonian coordinates, can be applied, so that the 
three spatial components of the geodesic equation become (Here we follow essentially the derivation in 29]) 

d 2 x i dx'dati dx> 1 

~w - - r ^^r " 2T ^ = ~ r44 = (2) 

where t denotes the Newtonian time. Therefore, we may apply a Newton's-like equations of motion for a scalar field 
<f> defined by 

dV = _ 90 

dt 2 dx i 1 ' 

Notice that <f> is not necessarily the Newtonian potential because the v « c condition was not applied to Einstein's 
equations. Comparing the above expression with (J2J, we obtain 

d(j) ldSh^ 

dx 1 ~ ~2 dx, ^ ' 

As the particle continues its free fall, while maintaining the slow motion, the gravitational field continuously builds 
up by small increments as 

g^v ~ r\^v + Sh^ + (Sh^) 2 H 

Actually, there is no reason to stop this process, so that (0J can be integrated for all perturbations of the Minkowski 
metric, up to a finite /i M „, leading to 



d(5h 4i ) = -_(i+ ff44 ) (5) 
z Jo 1 

This nearly Newtonian gravitational potential is characterized by an exact solution of Einstein's equations, with the 
non-linear effects present in the component (744 p9l |. 

In order to understand the implications of (JSj to the dark matter problem, suppose that we have the Schwarzschild's 
solution of Einstein's equations written in the usual spherical coordinates, so that 344 = — (1 — 2M/r). In this case, 
(JSJ) coincides with Newton's gravitational potential (f> — —M/r for a spherically symmetric gravitational source with 
mass M. If, this particular potential is applied to describe the motion of a star in a spiral galaxy corresponding 
to a spherically symmetric "visible mass" M, it does not describe correctly the rotation curves outside the galaxy 
nucleus, regardless of how strong the Schwarzschild field may be. On the other hand, if the star is close to the 
galaxy nucleus, then it will feel the pull of a spherically symmetric gravitational field which coincides with the above 
Newtonian potential —M/r. This coincidence explains why the two curves in Fig.l agree at the galaxy's nucleus. 
For any other solution of Einstein's equations which is not diffeomorphic to the Schwarzschild's solution, JSJ will 
produce a different near Newtonian dynamics. Our understanding is that the nearly Newtonian potential © carries 
a symmetry dependent non-linear effects contained in Einstein's equations through the component 544, as it will be 
exemplified in the next section. 

It is relevant to distinguish the present application of JSJ from a solution of the same problem using full general 
relativity as in |3Ct l3l|. Here, besides having lost general covariance as a consequence of the slow motion condition, 
only one component of metric has a direct contribution to the motion. In the following section we show that the 
velocity curves derived from (JSJ) for a vacuum gravitational field are compatible with the observed curves, using an 
exact solution of the vacuum Einstein's equations corresponding to a disk galaxy. 

III. VELOCITY CURVES NEAR A DISK GALAXY 

As a simple model for a disk galaxy we may consider a cylinder such that its height ho is much smaller than its 
radius tq. The line element produced by such object can be derived from the Weyl cylindrically symmetric metric, 
expressed in cylindrical coordinates (r, z, 9) as [33| 

dS 2 = e 2( - x -^dr 2 + r 2 e- 2 °d^ 2 + e 2 ^dz 2 - e 2 ° dt 2 (6) 
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where A = A(r, z) and a = a(r, z). The exterior gravitational field outside the cylinder, is given by vacuum Einstein's 
equations: 

-A, r + r<r 2 r - r<j 2 z = (7) 
— <j >r — ra \ rr — ra >zz = (8) 

\rr + \zz + <7 2 r + V% = (9) 

2ra tr a tZ = A jZ (10) 

To the above metric we apply the thin disk condition 

z e [-h /2,h /2], for r € [0,r ], h «r (11) 

In this case we may expand the functions a(r, z) and A(r, z) around z — 0, obtaining 

a(r, z) — a(r, 0) + za{r)... 
A(r, z) = \{r,0) + zb(r)... 



where we have denoted 

a{r) 



da{r, z) 
dz 



and b(r)= dX ^ z) 



d 



(12) 



z=0 



The thin disk condition l|ll|) implies that the above expansion can be truncated to the linear terms in z. Therefore, 
replacing a zz = and X :ZZ — in JSJl and i|l(J|) . they become simple ordinary differential equations on a, with general 
solution 

cr(r,z) = -T^lnr + c 2 (z) 

where we have denoted 

Ko = ^ (13) 
a(r) 

and where C2(z) is an r-integration constant. Derivation of a with respect to z gives C2(z) = a(r)z + Co, but since C2 
does not depend on r, it follows that a(r) — ao = constant. By similar arguments we find that b(r) = bg ^constant, 
so that Kq is also a constant. Replacing these results in Q and @, together with A )Z2 = 0, we also obtain ordinary 
equations for A(r, z). Therefore, the solution of Einstein's equations for the thin Weyl disk is 

a(r, z) — In r + a z + c (14) 
K 2 r 2 

Hr,z) = — Inr - a Q — + b z + d (15) 

where do is another integration constant. 

From (|14f) we obtain 1744 = — e 2<T = — e 2- ^ inr e 2a z e 2c _ Therefore, for an object moving in the disk plane z = 0, we 
obtain 

= -^(1 + .944)J z =o = -^(1 - e 2c V^) (16) 

In analogy with the Schwarzschild solution, the integration constant e 2c ° may be interpreted as a mass, with the 
difference that here we cannot compare it with the same Newtonian mass. However, we may assume that this constant 
is proportional to the visible mass M v of a disk-like galaxy (in units G=c=l): e 2c ° = 2/3qM v , where f3$ is a mass scale 
factor. It is even possible to interpret the factor /3 as something to do with the observed baryonic dark matter, but 
then we would require a correlation between the visible and dark matter in each galaxy |32| . 

The tangent velocity of a star as a function of the distance to the center of the galaxy can now be derived from 
the Newtonian- like equations of motion (jSJ), using the potential (J3J). Taking the force acting upon a star of unit mass 
with tangent velocity v — tur, ui — constant to be F — ^—f and comparing with the force generated by J^J F — —^r, 
we obtain, v 2 = so that for the considered disk we obtain ( in units G=c=l) 

u(r) = y/\PoM v K r K o\ (17) 
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The values of Kq given by (|13|) are determined by the coefficients of z in the thin Weyl disk metric. Interestingly, Kq 
is present even in the plane z — 0, a subtle consequence of the non-linearity of the vacuum Einstein's equations. 

In the thin disk case, the Newtonian velocity can be recovered for Kq = — 1 and [3q = 1. Since this particular value 
does not contribute to the rotation curves outside the galaxy core, the value Kq = — 1 must be ruled out for disk 
shaped galaxies. On the other hand, when Kq = +1, the velocity expression l|17f) does not correspond to any observed 
rotation curve. We conclude that \Kq\ must be smaller than one. 

Figure 2 shows the velocities calculated with (|17|l for a few known examples. Since the stars are supposed to be 
at the rim of each disk galaxy with radius ro, the origins of each curve were shifted, replacing r by r — ro, so that 
the shown curves start at the estimated ro for each galaxy. The values of Kq were determined by comparing the 
observed average velocity < vq > with the calculated velocity for each galaxy. In the given examples all values of Kq 
are positive but different, so that the curves have slightly different slopes. The values of (3q do not affect the shapes 
of the curves and were estimated for each galaxy from the know top speed in each case. 
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FIG. 2: Velocity curves from (1171 . for some examples. Distances are in Kpcs and velocities in Km/sec: 

(a) The continuous red line represents a simulation of the Newtonian curve for the Sun in the Milky Way. 

(b) Dotted red line is the Milky Way, for f3 = 1, K = 0.08, M v = 1 X 10 11 x M Q and r = 5 (at the Sun). 

(c) Magenta is NGC3198 for /3 = 1, K = 0.068, M v = 6 x 10 11 x M and r = 5. 

(d) Green is NGC3949 for (3 = 15.8, Ko = 0.13, M v = 2.5 x 10 9 x M and r = 1. 

(e) Blue is NGC3877 for f3 = 20, K = 0.18, M = 1.1 x 10 9 x M© and r = 3 

For comparison purposes we include below the observed plots (error bars) for NGC3877 and NGC3949 3] : 

UMa: NGC 3877 UMa: NGC 3949 
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FIG. 3: Observed rotation curves for NGC3877 and NGC3949 
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Summary 

The slow motion of objects in general relativity is described by the nearly Newtonian potential, obtained by imposing 
v << c in the geodesic equations only while leaving Einstein's equations and the geodesic deviation equations 
intact. The result is the nearly Newtonian gravity, something in between general relativity and Newtonian theory, 
characterized essentially by #44. The existence of such potential follows from fact that in general relativity the 
equations of motion are a consequence of the non-linearity of Einstein's equations, making a contrast with Newtonian 
gravity, where the equation of motion is postulated separately from the field equations. 

In particular, when the nearly Newtonian potential is derived from a vacuum solution of Einstein's equations, the 
slow motion of a test particle or a falling star is affected by the self interaction of the gravitational field, so that in 
principle there are no baryons involved. 

The loss of general covariance imposed by v << c means that the symmetry of the gravitational field solution of the 
vacuum Einstein's equations play a significant role in the velocity curves derived from (J5J, which is interpreted as a 
consequence of the non-linearity of Einstein's equations. In this respect, it should be noted that the Weyl cylindrical 
solution can be transformed to the Schwarzschild's solution by a diffeomorphism |34| . However, we cannot apply 
such transformation here because the diffeomorphism invariance has been lost. On the other hand, the solutions of 
Einstein's equations with a symmetry that resembles the gravitational field of a galaxy will describe velocity curves 
which are closer to the observed ones. This was exemplified by taking the Weyl solution with the format of a thin 
disk, as a model for a disk galaxy. In this case the velocity curves are remarkably close to the experimental curves. 

Clearly, a thin Weyl disk is a very poor mathematical model for a spiral galaxy. A more realistic model would be 
given by a static oblate spheroid, which can also be derived by a coordinate transformation of the Weyl metric |35| . 
Work on this is still in progress. 
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